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Abstract 

In this paper we review the asymptotic matching conjectures for r- 
regular bipartite graphs, and their connections in estimating the monomer- 
dimer entropies in d-dimensional integer lattice and Bethe lattices. We 
prove new rigorous upper and lower bounds for the monomer-dimer 
entropies, which support these conjectures. We describe a general 
construction of infinite families of r-regular tori graphs and give al- 
gorithms for computing the monomer-dimer entropy of density p, for 
any p G [0, 1], for these graphs. Finally we use tori graphs to test the 
asymptotic matching conjectures for certain infinite r-regular bipartite 
graphs. 
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1 Introduction 

The monomer-dimer covers of infinite graphs G, and in particular of the infinite 
graph induced by the lattice Z d , is one of the widely used models in statistical 
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physics. See for example [D U 12 U U O O UM CESl EH HH [20l [2H [23l [2H [251 

ESI ESI ED]. 



Let G = (V, E) be an undirected graph with vertices V and edges E. G 
can be a finite or infinite graph. A dimer is a domino occupying an edge 
e = (u, u) G -E". It can be viewed as two neighboring atoms occupying the 
vertices u, v G V and forging a bond between themselves. A monomer is 
an atom occupying a vertex w G V, which does not form a bond with any 
other vertex in V. A monomer- dimer cover of G is a subset £" of E such 
that any two distinct edges e, / G £" do not have a common vertex. Thus 
.E" describes all dimers in the corresponding monomer-dimer cover of G. All 
vertices V C V, which are not on any edge e G E', are the monomers of the 
monomer-dimer cover represented by E' . E' is referred to here as a matching. 
E' is called a perfect matching if V = 0, i.e. all the vertices of G are covered 
by the dimers. 

Consider first a finite graph G = (V, E). Then E' is called an l-matching 
if #E' = I. Note that 21 < #V. Let 0(7, G) > be the number of ^-matchings 
in G for any I G Z + . (Note that 0(0, GQ = 1 and 0(Z, (7) = if there are no 
/-matchings in G. Assume also that </>(/, G) = for a non-integer / > 0.) Then 
the monomer-dimer entropy of density p of G is defined as 



Let if)(x,G) := YH^Lo&ttityx 1 denote the matching generating polynomial of 
G. The pressure of G is defined as 



For an infinite graph G the monomer-dimer entropy of density p and the 
pressure Pcif) are defined by taking appropriate limsup on the finite sequences 
of graphs converging to G. (See for details §2.) 

We now consider the classical case in statistical physics: the lattice Z d , 
consisting of all d-dimensional vectors i = (ii, . . . ,i d ) with integer coordinates. 
(As usual we denote by Z,Z + ,N the set of integer, the set of nonnegative 
integers and the set of positive integers.) Let = (6ki, • ■ ■ , o~kd) be the unit 
vector in the direction of the coordinate Xk for k = 1, . . . , d. Then G(Z a! ) = 
(V = Z d , E), where (i, j) G E j - i = ±e fc for some k G [1 , d] . Note that 

G{lj d ) is an infinite 2d regular graph. 

Let h d (p) := h G{zd) (p) for any p G [0, 1] and h d := sup pg[01] h d (p). (h d and 

h d '■= h d (l) are called the <i-monomer-dimer entropy and the (i-dimer entropy 
respectively [12].) For d = 1 it is known that [121 §4]: 



Mp) = (i-f)log(i-|)-|log|-(l- P )log(l-p), pe[o,i]. (l.i) 



ha(p) 



logmax(0a^j,GQ,l) 



for any p G [0, 1]. 



Pait) : 



log^(e 2i ,G) 

w 
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The value of planar dimer entropy hz(l) was computed in [5] and [21] 



h 2 (l) = - > '—^ = 0.29156090 . . . 

K> 7r^(2g + l)2 

The exact values of ^(p) for p G (0, 1) and hd{p) for d > 3,p G (0, 1] are un- 
known. According to Jerrum |24j . the computation of the matching generating 
polynomials of finite planar graphs in general is computationally intractable. 
(This fact does not rule out the possibility that h d {p) are computationally 
tractable for d > 2, however for d > 3 it seems that hs(l) and /13 are hard to 
compute with high precision.) 

The properties of the entropy hd{p) for any p G [0, 1] was studied by Ham- 
mersley and his collaborators in [T5J [161 HZl US]- It was shown in [12] that hd{p) 
can be obtained from the limits of certain tori graphs, which are bipartite and 
2d regular. Using the proof of Tverberg's permanent conjecture, proved by 
the first name author [S], the following lower bound was shown in [T2] 



h d (p) > fh d (p) := |(-plogp-2(l-p)log(l-p)+plog2d-p) (1.2) 
for any p G [0, 1]. 

Tverberg's permanent conjecture states that the minimum of the sum of 
all Z x Z permanental minors of n x n doubly stochastic matrices is achieved 
only at the flat matrix J n — ( -) . It is a generalization of the van der Waerden 
permanent conjecture for doubly stochastic matrices, which is the case / = n. 
In [32] Schrijver gave a lower bound on the number of perfect matchings for 
r-regular bipartite graphs. It is an improvement of the lower bound implied by 
the van der Waerden permanent conjecture. Furthermore, this lower bound 
is asymptotically sharp. Equivalently, one can think that Schrijver's lower 
bound gives asymptotically the number of perfect matchings in large random 
r-regular bipartite graph. 

In [11] we stated a Lower Matching Conjecture, referred here as LMC, for 
any /-matchings of r-regular bipartite graph. For 21 = jfV this conjecture 
is asymptotically equivalent to Schrijver lower bound for perfect matchings. 
This lower bound can be viewed asymptotically as the number of /-matchings 
in a large random r-regular bipartite graph. The LMC implies the Lower 
Asymptotic Matching Conjecture stated in §2, referred here as LAMC, yields 
the following conjecture. 

h d {p) > gh 2d (p), for any p G (0, 1] and d > 2, (1.3) 

where 

1 J) 

gh r {p) := - (plogr -plogp- 2(1 - p) log(l -p) + (r -p) log(l - -)), (1.4) 

for any integer r > 2. Note that hi(p) = gh,2(p). In a recent paper [10] 
the LAMC was proven for the sequence of densities p = s = 0, 1, . . ., for 
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any given r > 2. Hence (11 .3p holds for p = ^frj , s = 0, 1, . . .. In particular 
Ml) > 0^(1) for any d G N. The inequality ft 3 (l) > 0.440075 is the best 
known lower bound. A recent massive computation performed by the third 
named author in [26] gives the best known upper bound hs(l) < 0.457547. 

The conjectured lower bound (11.31) yields a lower bound for the <i-monomer- 
dimer entropy h d . In particular, the conjectured lower bound (I1.3P yields 
h 3 > 0.784992989. The validity of QEJD for d = 3,p = ^ = § implies the 
best lower bound known /i 3 > /13(f) > .7845241927. In this paper we give 
new lower bounds on hd(p) which yield the inequality h 3 > /i 3 (0.6814) > 
.7849602275. The numerical computations in [T2] yield the best known upper 
bound h 3 < 0.7862023450. 

In [TTJ we stated an Upper Matching Conjecture, referred here as UMC. 
Namely, let K r>r be a complete bipartite graph on 2r vertices, where the de- 
gree of each vertex is r. Denote by qK rr be the graph consisting of q copies 
of K r r . Then the UMC claims that any r- regular bipartite graph G on 2qr 
vertices satisfies <fi(l, G) < <p(l, qK r ^ T ) for I = 0, . . . , qr. We also have a cor- 
responding Upper Asymptotic Matching Conjecture, referred here as UAMC, 
which is slightly more technical to state. (See §6.) For r = 2 we proved these 
conjectures in (TTJ. 

The main purpose of this paper is to give theoretical and numerical evi- 
dences on the LAMC and UAMC and their applications to the estimates of 
the monomer-dimer p-densities for Z d and for the Bethe lattices, i.e. (i-regular 
infinite trees. We believe that the computational and theoretical setting dis- 
cussed in this paper are of interest by itself and to researchers in asymptotic 
combinatorics, which is widely used in statistical physics. 

We now outline briefly the main setting of our computations for the verifi- 
cation of the two asymptotic conjectures. It is well known that the asymptotic 
growth of many configurations in statistical physics are given in terms of the 
spectral radius of the transfer matrix. See for example [12]. In this paper we 
construct infinite families G n = (V n , E n ),n G N of r-regular bipartite graphs, 
which are coded by a specific incidence matrix A G {0, \y NxN . This sequence 
of graphs converges to an infinite r-regular graph G. Using programs based on 
software developed by the third named author one obtains the transfer matrix 
B(t) G M 2 x2 , corresponding to the matching generating polynomial with 
the value x = e 2t . Since the infinite tori graphs corresponds to a subshifts of 
finite type, abbreviated here as SOFT, one can compute the pressure function 
P(t) in terms of the spectral radius p(B(t)). This is well known to the ex- 
perts, and we bring the proofs of these formulas in the paper for completeness, 
using the general techniques in [13] . (The properties of the pressure function 
P(t),t G K, for mult i- dimensional SOFT, as for example the monomer-dimer 
models in Z d ,d > 1, are studied in detail in [T3]-) Then the monomer-dimer 
]?-density hc(p) is computed by using p(B(t)) and its derivative. (In this set- 
ting p = p(t).) We then compare ho{p{t)) to the upper and lower bound given 
by the lower and upper asymptotic conjecture. 
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We now briefly survey the contents of our paper. In §2 we discuss the 
monomer- dimer entropy h{ Gn y(p) of density p and the pressure function P{G n }if) 
for a sequence of finite graphs {G n = (V n , E n )} of bounded degrees such that 
jfV n — > oo. We define the function low r (p),p G [0,1] which gives the sharp 
inequality /i{G„}(p) > low r (p) for any sequence of bipartite r-regular graphs 
and any p G [0, 1]. We state the LAMC, which is equivalent to the equality 
low r = gh r . Furthermore if the sequence {G n } is a sequence of random r- 
regular bipartite graphs we conjecture that hiQ\(p) = gh r {p) almost surely 
[TT] . In §3 we use the recent verification of the LAMC for p = s = 0, 1, . . . 
for any r > 2 to derive tight lower bounds on for low r . In §4 we discuss the 
applications of our results to Bethe lattices, i.e. infinite dimensional r-regular 
trees. In §5 we discuss the sequence of tori graphs, which are considered in 
[T2] and [13] to compute h 2 , h 3 and h 2 (p). We prove the thermodynamics for- 
malisms for such graphs which gives the monomer-dimer entropy of density 
p in terms of the pressure. In §6 we describe a fairly general construction of 
sequences of regular graphs, which includes the sequence of tori graphs. In §7 
we describe the upper matching conjecture and its asymptotic version, called 
the upper asymptotic matching conjecture. We give upper bounds for h{a n }ip) 
for any sequence of bipartite r-regular graphs and show that in some regions 
these bounds are relatively close to the UAMC. In §8 we describe our com- 
putational results, which support the conjectures stated in this paper. In §9 
we identify an infinite graph with the maximal pressure among other infinite 
graphs in certain families of sequences described in §6. 



2 Entropies, Pressure and LAMC 

We will now define a limiting monomer-dimer density for a sequnce of bounded 
degree graphs. 

Definition 2.1 Let {G n },G n = (V n ,E n ),n G N be a sequence of finite 
graphs, where multi edges are allowed, such that jfV n — > oo and the degree of 
each vertex in G n is bounded by d for n G N . Forp G [0, 1] we define h{G„}{p), 
the monomer-dimer entropy of density p, as follows: 

7 /x r log (p{l n ,G n ) 
h{G n }[P) = hmsup — , (2.1) 

n— »oo if Vn 

21 n 

over all sequences l n G Z + satisfying lim — — =pg [0,1]. (2.2) 

h{G n }iX) an d ^{G„} := su P P g[o,i] h{G n }(p) are called the dimer entropy of {G n }, 
and the monomer-dimer entropy of {G n } respectively. For t G H. the pressure 
of {G n } is defined as 

P{G n }{t) := hmsup — . (2.3) 
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Let G = (V, E) be an infinite graph, where multi edges are allowed. Assume 
that the maximal degree of vertices in G is d < oo. A sequence of graphs 
{G n }, G n = (V n , E n ), V n C V, n = 1, 2, . . . , V n , where multi edges are allowed, 
converges to G if the following conditions hold: 

1. V\ C V<2 C . . . are finite subsets of V satisfying the condition V = 

u~ xK = v. 

2. Each G n contains the induced subgraph G(V n ) := (V n ,E(V n )) of G on 
the set of vertices V n , and the degree of each vertex in G n is at most d. 

3. Let v G V n and assume that all neighbors of v in V are in V n . Then E 
and E n have the same set of edges that contain v. 

Then h G (p) := h {Gn} (p), h G := h {Gn} , P G (t) := P{ Gn }(t). 

The above definition of entropy and pressure of an infinite graph G depends 
on the specific choice of the convergent sequence {G n } to the infinite graph G. 
For G(Z d ) one has a whole class of the sequences {G n }, for which the resulting 
h{G n }(.P)i P{G n }(t) is independent of the choice of the convergent sequence {G n } 
[T5l fT6l fT71 fT9l HBj. In this case we denote by hd(p), Pd(t) the corresponding 
quantities. For other infinite graphs G discussed in this paper we choose a 
convenient convergent sequence {G n }, and we do not discuss the corresponding 
class of sequences which yield the same entropy and pressure. 

The properties of the entropy h^p) for any p G [0, 1] was studied by 
Hammersley and his collaborators in [15] [TBI E33 EE]- Let us mention two 
properties that are of interest in this context. For any m 6 N let (m) := 
{1, . . . , m} = [l.m] n Z be the set of integers between 1 and m. For any 
m = (mi, ■ ■ ■ , fUd) € N d let (m) := (mi) x . . . x (rrid) C N d is the set of points 
in the lattice Z d located in the box [1, mi] x . . . x [l,md] in M. d . Denote by 
vol(m) := Ylt=i m i the volume of the box (m). Let G(m) := ((m) , E(m)) be 
the subgraph of G(Z d ) induced by (m), i.e. (l, j) G E(m) <^=^> l, j G (m) and 
i — j = ±ek for some k G (d). Let m n := (mi jn , . . . , md, n ), n G N be a sequence 
of lattice points in N d , such that m n — > oo <^=^ — > oo as n — > oo for 
each k = 1 , . . . , d. Then for any sequence /„ G [0, ^p^] n N the following 
conditions hold: 

log0(/ n ,G(m n )) 2/ n 
hm = h d (p) if m n —> oo and hm — = p G 0, 1 . 

n^oo VOl(m n ) n^oo VOl(m„) 

(2.4) 

The above characterization yields that hd(p) is a concave continuous function 
on [0,1], see p]. 

Let T(m) := ((m), ^(m)) be the torus on (m). Thus two vertices i, j G (m) 
in T(m) are neighbors if (i,j) G E(m), or for any m^ > 2 the vertices 
(ii, ... ... ,^) and (i x , . . . , m k , ik+i, ■ ■ ■ , id) are adjacent for 

any k G (d) and (i u . . . , z fe _i, i k+1 , . . . , i d ) G N^ 1 . Clearly 0(/,T(m)) > 
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4>(l,G(m)) for any / £ Z + . It was shown in [12] that the condition (12.41) 
can be replaced by the corresponding condition on the torus: 



lim 



log</>(/ n ,T(m n )) 



h d {p) if m 



oo and lim 



21 



ri 



P e [0,1]. 

(2.5) 



vol(m n ) 



vol(m n ) 



n 



oo 



(It is assumed that l n £ [0, ( ^" J ] n N.) More general, one can show that 
hd(p) = h G(zd) (p). 

There are several advantages of considering T(m) over G(m). Assume 
that rrik > 2 for k = 1, . . . , d. First, the graph T(m) is a 2c/-regular graph. 
Second, the automorphism group of T(m) is quite big, which can be very well 
exploited, using the general method of [26]. See also [27]. and [121 13] f° r the 
computations of h d , h d and h d {p) respectively. 

The fact that T(2m) is 2o?-regular bipartite graph was exploited in [12] 
to show (11.2)1 . This lower bound is obtained by noting that if G = (V,E) 
is an r-regular bipartite graph then <p(l,G) > f r (l,#V), where the function 
f r (l,2n) := (?) W (n) is determined from the proof of Tverberg's permanent 
conjecture [8]. 



The LMC stated in [11] claims that <f>(l, G) > g r {l,#V) (> f r (l,#V)) for 



For 21 = jfV (I2.6P is Schrijver's lower bound for perfect matchings in r-regular 
bipartite graphs on 2n vertices. The LAMC, which yields fll.3[) . can be stated 
as follows: 

Conjecture 2.2 ( The Lower Asymptotic Matching Conjecture.) 

Let Q(2n, r) be the set of r-regular bipartite graphs on 2n vertices, possibly with 
multi edges. For each I £ [0,n] PlZ let a(l,2n,r) := minG e g( 2n ,r) <p{l, G). For 
p £ [0, 1] let \ow r (p) be the infimum liminf^oo log At( -^ rifc ' r - > over all sequences 
< Ik < Tiki k £ N such that lirrik^ooT*' = V- Then 



The results in [TT] show that for a given p £ [0, 1] and r > 2, the above 
conjecture is equivalent to the statement that the number of /-matching in a 
random bipartite r-regular graph will behave asymptotically as in Conjecture 
12.21 In particular, the random graphs minimize, in the asymptotical sense, the 
number of /-matchings in r-regular bipartite graphs. 

It is shown in [TT] that the LMC and LAMC hold for r = 2. Furthermore, 
the cycle C 2n on 2n vertices satisfies the inequality </>(/, C 2n ) < </>(/, G) for any 
G £ £(2n,2). Hence 



any r-regular bipartite graph, where. 




(2.6) 



low r (p) = gh r {p). 



(2.7) 



k— »oo 



lim 



log0(/ fc ,C 2TOfc ) 
2n k 



gh 2 (p) = hip) 
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for any sequence < l k < n k ,k E N satisfying the assumptions of Conjecture 

In a recent paper [TUl Theorem 5.6] the following results were proven: 

Theorem 2.3 The Lower Asymptotic Matching Conjecture holds for the 
following corresponding sequence of densities p = = 0,1,..., for any 

given r. In particular lil.3\) holds for p = gjzjj s = 0, 1, . . .. 

This can be extended to give a bound for all p in the following way. 

Definition 2.4 For 2 < r E N let ghl r (p),p E [0,1] be the following 
function. 

• 9hl r (^- s ) = gh(^r s ) for s = 0, 1, . . .. 

• ghl r (p) is linear on the interval [ / , x , for s — 0, 1, 

• ghl r (0) = 0. 

The concavity of hd(p) and Theorem 12.31 yields: 

h d {p) > ghl 2d {p), for all p E [0, 1] and d = 2, 3, ... . (2.8) 

In the next section we improve substantially these lower bounds. 

In [TOl Figure 1] are plotted the graph of gh^p), the graph corresponding 
to UAMC and the 19 values of the /i2(p) computed by Baxter pQ. (Baxter's 
computations are based on sophisticated heuristical arguments. His compu- 
tations were recently verified by rigorous mathematical methods in [13].) It 
turns out that Baxter's values are very close to the values of gh±(p). 

In Figure [1] we show the graphs of ghl r (p), low ri i(p), a lower bound for 
low r (p) given in the next section, and gh r (p) for r = 4. Note that the dif- 
ferences of the three graphs are relatively large on the first interval from the 
right [^q-j-, 1], slightly less on the second interval from the right [^, ^-], and 
ignorabie from the fourth interval to the right [^rg, ^g]- We notice that the 
differences between the functions ghl r (p), low ri i(p), gh r {p) decrease as r in- 
creases. (This observation applies also for the values r = 3, 6 which are not 
plotted here.) 

3 Lower bounds for low r (p) 

In this section we give a lower bound for the function low r (p), which is defined 
in Conjecture 12.21 

Theorem 3.1 The function low r (p) + |(plogjo + (1 — p) log(l — p)) is 

concave. 
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Proof. Let G G Q(2n,r). Consider the polynomial 9(x) := (—x) n ip(—-,G). 
Since 0(7, G) > for I = 0, . . . , n it follows that has n complex nonzero 
roots. It is well known [2D] that 9(x) has only positive roots. The Newton 
inequalities, see e.g [29] . yield 

4,(1 -l,G)4,(l + l, G) , 4>(l,G) ,2 . 

Let Gi t 2n,r £ (?(2n, r) be an r-regular graph for which the equality //(/, 2n, r) = 
4>(l,Gi } 2n,r) holds. (13.11) and the minimal characterizations of fi(k,2n,r),r = 
0, . . . , n yields 

H(l-l,2n,r)fj,(l + l,2n,r) f fi(l,2n,r)) ^ 2 , , * 
G5 © = (3-2) 

This is equivalent to the statement that the sequence 

Oi,2n,r : = log ^"4^" ^ °> i = 0, . . . , n 

is a concave sequence. Let 2n, r) be a piecewise linear function on [0, 1] 
defined as follows: 

. a (I,2n,r) = 2f* Z = 0,...,n. 

• a(-, 2n, r) is linear function on the interval [-, — 1 for / = 0, . . . , n — 1. 

The concavity of the sequence ai^ n ,r, / = 0, . . . , n is equivalent to the concavity 
of a(x, 2n, r). Let {nk}^? =l be an increasing sequence of positive integers. Let 
h € [0,rifc] fl Z, fc G N be a sequence satisfying lim^oo ^- = p 6 [0,1]. Use 
Stirling's formula to deduce that 

log i 

lim fc = --(plogp+ (l-p)log(l -p)). (3.3) 

fc^oo zrik z 

It is straightforward to show that 

liminf a(x, 2n, r) = low r (x) + - (plogp + (1 — p) log(l — p)) . (3.4) 

Since each a(-, 2n, r) is concave it follows that \ow r {p) + \ (plogp+(l— p) log(l — 
p)) is concave. 

□ 

The arguments of the proof of the above Theorem combined with the def- 
inition of hd(p) = h G (id)(p), implies a stronger concavity result than given in 

m. 
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Corollary 3.2 Let hd(p),p G [0, 1] be the monomer- dimer entropy of den- 
sity p for the graph G(Z d ). Then h d (p) + ^(plogp + (1 — p) log(l — p)) is 
concave. 

Theorems 12.31 and 13.11 yields: 

Corollary 3.3 Let low rj i (p),p G [0, 1] be defined as follows. 

• low r ,i(^) = 9 h r(^- s ) for s = 0, 1, . . .. 

• low rl (p) + 7|(plogp+(l — p) log(l— p)) is linear on the interval [ r , r s+1 , 
for s = 0, 1,. . .. 

• low r>1 (0) = 0. 

Then low r (p) > low rj i(p) for any p G [0, 1]. 

Figure [1] shows the position of the graphs ghl r (p) < low rj i(p) < gh r (p) for 
r = 4. 

We now give a different lower bound for low r (p) using [10j Theorem 5.6]. 

Theorem 3.4 Let low r ^{p),p G [0, 1] be defined as follows. 

• low ri2 (^) = gK(^) for s = 0, 1, . . .. 

• For p G (^tj, 1) low rj2 (p) the maximum between the two following num- 
bers 

p 11 

-(logr + (r - l)log(l - -)) - -(plogp + (1 -p)log(l - p)) 
and 

p r 1 

■ logr — plogp — (1 — p) log(l — p) H — log(l — 



2 ° r OL y r/ ov r ' 2 ov r + r 

For p G ( r+ g +1 ; 7^7) low rj 2(p) i/ie maximum between the two following 
numbers 

p 1 

-logr + - ( -jo log jo - 2(1 - jo) log(l - jo)) + 

\ ( (r + s - 1) log(l - — j— ) - (s - 1 + p) log(l - ] —^\ 
2 \ r + s s 

and 

p 1 

- logr + - (-plogp - 2(1 - p) log(l - p)) + 

1 / 1 1 -p 

- (r + s) log(l — ) - (s + p) log(l — 

2 \ r + s + 1 s + 1 

for s = 1, . . .. 
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• low rj2 (0) = 0. 

Then low r (p) > low rj2 (p) for any p G [0, 1]. 
Proof. [TDJ Theorem 5.6] states 

liminf °^ ^ kl ^ fc ) > P \ Q g r _| — (—p\ogp — 2(1 — p) log(l — p)) + (3.5) 
n— »oc 2rifc 2 2 

i ( (r + j - 1) log(l — ) -(j-l+p) log(l ' ~ ; '' 



2 V r - j ' j 



for any sequence of Gf~ G Q(2n,k,r) where lim^oo ^ = p and any j G N. 
This implies the the inequality low r (p) > low rj2 (p) where p G [ T . + ^ +1 , ^rrj] and 
s G N. (Choose j = s, s + 1.) 

In the interval [^fj, 1] the second inequality follows from the above inequal- 
ity for j — 1. For the first inequality we use the arguments of Theorem 13.11 
Combine the arithmetic-geometric inequality and Schrijver's inequality [32] to 
deduce 

i i (r-1)- 1 



r-2 



(?) ' 

Taking logarithm of both sides, dividing by 2n and using any sequence 
satisfying (I2.2p we deduce the inequality 

p 11 
low r (p) > -(logr+(r-l)log(l--)) - -(plogp+ (1 -p) log(l -p)) (3.6) 

for all p G [0,1]. □ 

It turns out that for many of the values of p, the lower bound \ow ri2 (p) a 
better lower bound than low n i(p), and it is very close to the function gh r (p). 
Figure [2] compares the differences low n i(p) — gh r (p), plotted in black, and 
low n2 (p) — gh r {p), plotted in blue, for r = 4. 

From this graph and the graphs for r = 3,6, which are not plotted here, we 
conclude that the errors gh r (p) — ghl r (p), gh r (p) — low n i(p), gh r (p) — low ri2 (p) 
are decreasing monotonically with r. 

Since k>W6, 2 (0.6814) = .7849602275 we deduce the inequality 

h 3 > /i 3 (0.6814) > low 6 (0.6814) > low 6 , 2 (0.6814) = .7849602275 
given in §1. Combine Corollary 13.31 with Theorem 13.41 to deduce 
Corollary 3.5 For 3 < r G N and p G [0, 1] 

low r (p) > max(low r x(p), low r 2 (p)). 
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4 Monomer- dimer densities for Bethe lattices 



Let T(r) be an infinite r-regular tree. Recall that T(3) is known the Bethe 
lattice. Clearly, each T(r) is bipartite. For each T(r),r > 2 we construct a 
convergent sequence {G n (r)} in sense of Definition 12. 11 

Fix a vertex O in T(r) and consider all vertices in T(r) whose distance from 
O is n > 1. Then the number of such vertices is r(r — l)* 1-1 . Let Ai> . . . , A r 
be the r vertices of distance 1 from 0. The number of vertices in T(r) whose 
distance from O is exactly n is divided to r classes Ai yU) i = 1, . . . ,r, where 
the points in Ai )U have distance n — 1 from Aj. Let X n = U^Ai^. Note 



that #X n = #A r , n+1 = (r - 1)™. Let 7„ := {0} U^ = i A,j U 4-,n+i Let 



H n (r) = (X n U Ar^+i, F n ) be an arbitrary r — 1 regular bipartite graph with 
the two classes of vertices X n) A r>n +i. Let G n {r) = (V n ,E n ) where E n are the 
union of the edge set in the induced graph T(r)(V n ) and the set F. Note that 
G n (r) is r-regular and bipartite. Then G n (r), n = 1, . . . , converges to T(r). 

Note that T(2) is isomorphic to the integer lattice Z, and G n (2) is a cycle 
of length 2{n + 1) for n G N. Hence 



In [TT] it is shown that k>W2(p) = gh 2 (p). 

Using the definition low r (p) and the fact that G n {r) are r-regular and 
bipartite we obtain. 

Corollary 4.1 Let 2 < r G N. TTien 



For more complex models, like the Ising model, it is known that this kind 
of limit is sensitive to the exact limiting sequence of graphs [18] . It is an 
interesting problem if equality holds in the above inequality for some choices 
of random graphs H n (r) G <?(2(r — 1)™, r — 1), n G N. 

5 An example of sequence of tori 

We first discuss a sequence of graphs that give the lower and upper bounds 
for hd and hd(l) for the graph G(Z d ) considered in [T2J. Assume that the 
dimension d > 1. Let m' := (mi, . . . , md-i) G N d_1 be fixed and assume 
that rrii > 3 for i = 1, . . . , d — 1. Consider the sequence of ci-dimensional 
tori T((m',n)) = (V n ,E n ),n = 3,4, . . .. Each torus is a 2d-regular graph. If 
mi, . . . , rrid-i and n are even then T((m',n)) is bipartite. The vertex set of 
T((m / , n)) is the set V n := ((m', n)). V n can be viewed as composed of n layers 
of vertices (m'). The edges between all vertices (m') in each level k are given 
as in the d — 1 dimensional torus T(m'). The other edges of T((m',n)) are 
going from level k to level k + 1 for k = 1, . . . , n, where the level n + 1 is 
identified with the level 1. (We also identify level with the level n.) The rule 



ha(Z){p) = h {Gn{2 )}(p) = gh 2 (p) for all p G [0, 1]. 



(4.1) 



h{G n (r)}(p) > low r (p) for allpe [0, 1]. 



(4.2) 
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for the edges between the level k and the level k + 1 is independent of k. Thus 
the vertices (Y, k) and (j', k + 1) in V n are adjacent if and only if i' = j'. The 
adjacency matrix between the two vertices i' in the level k and j' in the level 
k + 1 is given by the — 1 matrix A{ra!) := (ai'j')i'j'e(m')j which is an identity 
matrix of order iV = vol(m'). For any square matrix A G R nxn we denote by 
tr A and p(A) the trace and the spectral radius of A respectively. 

Let us recall first the computation of the monomer-dimer entropy hd given 
in [12]. The entries of the transfer matrix B(m') = {bsT)s,Tc{m') are indexed 
by two subsets of S, T of (m'). (These subsets may be empty.) First b ST = 
if SnT 7^ 0. Second assume that Sf)T = then bsr counts the number of the 
monomer-dimer covers of the subgraph of T{va!) induced by the set vertices 
(m')\(5 U T). Note that any subgraph of T(m') induced by a set U C (m') 
can be covered by monomers. Hence bsr > 1- (If 5 U T = (m) then b$T — 1-) 
It is not hard to see that the product of n terms bs 1 s 2 bs 2 s 3 ■ ■ ■ &s n _iS„&s n Si 
corresponds to all monomer-dimer covers of G n = T((m', n)) with the following 
conditions. For each level k — 1, . . . , n the dimers going from the level k to 
k — 1 are located at the set Sk and the dimers going from the level k to the 
level k + 1 are located at the set Sk+i- Let $(G> t ) := YlHo <K^> G n ) number of 
all possible monomer-dimer covering T((m',n)). Then tr B(m') n = $(G n ). It 
is shown in |12| 



Inn '°^' G -» = ^0[ B W), hd = lim (5.1) 

n^oo jfV n Vol(m') m'^oo vol(m') 

logp(£(2m')) , ^ d ! 

and h d < \ V , ;; for any m' G N . 
vol(2m'J 

(Here 2m' := (2m l5 . . . , 2m,j_i).) The lower bounds for h d are also expressed in 
terms of linear combinations of certain log p(B(m')) corresponding to different 
values of m'. 

Let T((m',Z)) be an infinite graph given by the set of vertices (m') x Z 
and the following set of edges E(m',Z). ((i',p), (j', q)) G E(m',Z) if either 
p = q and (i',j') G ^(m') or \p — q\ = 1 and i' = j'. Thus the sequence 
of graphs G n := T((m',n)),n = 3,4,... converges to G := T((m',Z)). Let 
ho(p) '■= h{G n }(p) be defined by (j2.lH2.2p . We now show how to compute 
ha(p) using the pressure function. 

Let S,T be two disjoint subsets of (m'). Let E' C ^(m') be an /-matching 
of T(m') so that each edge (u, v) G E' represents a dimer occupying two adja- 
cent vertices in T(m') located in (m')\(SUT). To this matching we correspond 
a monomial x . Let cst{%) be the sum of all such monomials. cst{x) is the 
matching polynomial for the graph T(m',S, T), the subgraph of T(m') in- 
duced by the subset of vertices (m')\(5UT). We let c ST (x) = if S, T C (m') 
and S n T ^ 0. Then b ST = c ST (l). Let 6 5 t(^) := c ST (e 2t )e^ s+ # T ^ and 
B(m, t) := (bsr(t))s,Tc{Ta>)- The arguments in [T2J that show that $(G n ) = 
^(l,G n ) yield the equality tr5(m,t) n = ^(e 2 *,G„). The definition (TOD of 
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pressure P(t) and the arguments in [12] for the equality (15.11) imply 

PM;= q c . )W = ^K')) , ieR . (5.2) 

(We suppressed the dependence of P(t) on m'.) 

The following results are known, e.g. [U [T3], and we bring their proof for 
completeness. 

Theorem 5.1 Let P(t) := P{g„} be defined by Ii5.2\) . Then P is a smooth 
increasing convex function on R. Furthermore 

lim P'(t) = 0, lim P'(t) = 1. (5.3) 

Let p{t) := P'(t),t G R. TTjen p(£) zs an increasing function such that 
p([— oo,oo]) = [0,1]. Furthermore 

h{G n }(P'(t)) = P(t) - tP'(t) for any t E R, (5.4) 

and h{c }{p) is a continuous concave function on [0,1], which is smooth in 
(0,1). 

Proof. The well known result [22] yields P(t) is a convex function of 
t. Since m' is fixed we let B(t) := B(m',t) and iV = vol(m'). As B(t) is 
irreducible, and the nonzero entries of B(t) are increasing function on M, it 
follows that P(t) increases. Since p(B(t)) is a simple root positive of the char- 
acteristic polynomial it follows that p(B(t)) and P(t) is an analytic function 
in some open domain containing K. 

N 

B(t) = ^^e u Bi, each Bi is nonnegative. (5-5) 

Bo = (&sT,o)s,Tc(m') corresponds to monomer configurations of . I.e. B Q is a 
diagonal matrix with one nonzero entry = 1. The matrix B^ corresponds 
to the tiling of T((m',n)) by dimers. Hence P(t) = log p(B(t)) is not a 
linear function. The analyticity of P(t) yields that P"(t) may have only finite 
number of zeros on any closed interval [a, b}. The convexity of P(t) implies 
that P" > on R. Hence P" is positive on any [a, b] except a finite number 
of points. Thus p(t) = P'(t) increases on R and Pit) is strictly convex on 
R. Let x = e*. Then the B(x) := B([ogx) is a polynomial in x. Since 
p(B(0)) = p(Bq) = 1 is a simple root of the characteristic polynomial of B it 
follows that p(B(x)) is analytic in some disk |x| < e, such that \p(B(x))\ > 
in this disk. Hence the branch log p(B (x)), log p(B(0)) = is analytic in this 
disk and has Taylor expansion. The same statement holds for the derivative 
of log p(B(x)). Substitute x = e* to deduce that P(t) = logp ^ e u anc l its 
derivative have convergent series in x = e* for t < —R, for some R » 1. This 
implies the first equality in (15.31) . 
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Observe that 

Hence P'(t) = 1 + ( logp( ^^ )' . The arguments above for the first equality in 



( 15. 3p imply that the second equality in (15. 3p . 
We now show the inequality 

P(t) >pt + h G (p) for all p G [0, 1]. (5.6) 

Let l n G [0, fi]nZ,n6N be sequence satisfying ( 12. 2ft . Then 

ti B(t) n = ^e 2t ,G n ) > e^^Gn) => 
log tr B(t) n > 2U + log0(/ n ,G ri ) 



Nn ~ Nn Nn 
Recall that #V n , the number of vertices in G n is Nn. Use the definition of 
P(t) and (T22D to deduce 

P{t) > pt + hm sup — — . 

Use the definition ( 12. 1112. 2p of h G (p) '■= ^{G„}(p) to deduce the inequality ( 15. 6L 
It is straightforward to show that h G (p) upper semicontinuous on [0, 1]. 
We now show that for each t G 1R there exists p(t) G [0, 1] such that: 

P(t)<p(t)t + h G (p(t)). (5.7) 

Let l n (t) G [0, ?f] n Z] satisfy 



■in 



nN 



Hence 



< ™<l>{l,G n ) < e a »WVa«(*),Gf n ), for / = 0, . . . [— J. 

</>(e 2t , G„) < ye^^^at),^). (5.8) 

Take a subsequence n^, k G N such that lirn^oo log ^ e ' G " fc ^ = P(t). Choose 

a subsequence jfc, fc 6 N of n^, A; G N, such that lini/^oo -^-^ = p(t) G [0, 1]. 
Take the logarithm of the inequality (15.81) and divide by nN. Let n = jk and 
let k — > oo. The definition (12.1112.21) of h G (p) yields the inequality (15. 7p . 

The inequalities (15.61) and (15.71) yield the equality P{t) = p{t)t + h G {p{t)). 
Moreover (15.61) yields P(y) lies above the line p{t)y + h G (p(t), which intersect 
P(y) at the point y = t. Hence p(t) = P'(t) and h G (p(t)) = P(t) - p(t)t. 
I.e. ( 15.41) holds. Since P' increasing and analytic the implicit function the- 
orem yields that t = Q{p) is analytic in p G (0,1). Hence h G {p) is analytic 
on (0, 1). Observe that —h G (p) is the Legendre function corresponding to a 
smooth strictly convex function P(t) [31]. Hence h G (p) is concave on (0,1). 
Our arguments yield that h G (p) is continuous on [0, 1]. Hence h G (p) is a con- 
cave function on [0, 1]. □ 
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Remark 5.2 Let hc{p) '■— ^{g„}(p) an d P(t) := P{G„}{t) be given as in 
Definition \2. 11 Theorem \5.1\ applies to hc{p) and P(t) in the following cases: 



1. There exists a nonnegative irreducible matrix B(t) of the form ( 15.51) such 
that 

• #V n = nN, ne N. 

• ^(e 2 *,G„) = tvB(t) n , n G N. 

• p(B> ) = 1 and p{B^) are positive simple roots the characteristic 
polynomials of B and B^ respectively. 

2. G n is a disjoint union of n copies of a finite graph H = (W, F) which 
has a perfect matching. Then P(t) = log t^ — ■ 



6 A construction of sequences of graphs 

We now generalize the construction in the previous section to a general con- 
struction of a sequence of regular graphs. Let F — (U, D) be an undirected 
graph with the set of vertices U and the set of edges D. For n > 2 let 
Gn '■= (V n , E n ) be the following graph. V n = U x (n), i.e. we can view V n 
consisting of n copies of U arranged in the n layers (U, 1), (U, 2), . . . , (U, n). 
We let {U, 0) := {U, n), (U, n + 1) := (U, 1). Then 

1. For any u,v G U and k G (n) ((u, k), (v, k)) G E n <^=^ (u, v) G D. 

2. Any other edges of E n are between the vertices (U, k) and (U, k + 1) for 
k = 1, . . . ,n. 

3. Let A = (a uv ) UtVe u be a given nonzero — 1 matrix. Then for each 
k G (n) ((«, k), (v, k + 1)) G i?„ <^=^> a uv = 1. We call A the connection 
matrix. 

4. For any two subsets S,T C U, (S,T may be empty), let cist G Z + be 
defined as follows. If #S ^ #T then clst = 0. Assume that #S = #T. 
Let i^S 1 , T) be the set of all bijections j3 : S T. Then = 1 and 

&ST = J2/3eB(S,T) UsGS a si3(s) for #S = #T>1. 

Thus a^T is the number of perfect matchings in the subgraph of the 
bipartite graph on the set of vertices (U, 1)U(U, 2), and the set edges E C 
(U, 1) x (U, 2) given by A, and induced by the subset of vertices (5*, 1) U 
(T, 2). Let A := {ost)s,tcu be a 2* 17 x 2* 17 matrix with nonnegative 
integer entries. 

5. For any two disjoint subsets S,T C U, let cst(x) be the matching gen- 
erating polynomial of the subgraph of F induced by the set of vertices 
U\(S U T). For non-disjoint subset S, T C U let cst(%) = 0. Let 
M(t) := (c ST (e 2t )e(# s+ * T ») s ,Tcu and S(t) := M{t)A be 2* 17 x 2* 17 
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nonnegative matrices for any i e K. Then log p(B(t)) is a continuous 
convex function on IR. If B(l) is an irreducible matrix then log p(B(t)) 
is an analytic function on 1L (See arguments of the proof of Theorem 



Then the sequence G n , n = 2, . . . has the following properties: 

• If F is connected then each G n is connected. 

• Assume that F is bipartite, where U = U\ U U 2 , D C U\ x U 2 . Suppose 
that the edges between the two consecutive levels of vertices (U, k) and 
(U, k + 1) are either between (Ui, k) and (Ui, k+ 1) for i — 1, 2 or between 
(Ui, k) and (U i+ i, k + 1) for i — 1, 2. (f/ 3 := L^.) If n is even then G n is 
bipartite. 

• Assume that F is jo-regular. Assume that the matrix A has q l's in each 
row and column. Then G n is p + 2g-regular graph. 

• Assume that F is p-regular bipartite. Let U — U\ U £/ 2 ,-D C f/i x U 2 
and n is even. Assume that the matrix A has the following properties. 
Each row indexed by u G U\ and each column indexed by v G U 2 has q 
l's, and each row indexed by v G £72 and each column indexed by u G £7i 
has g — 1 l's. (q G N.) Then G n is p + 2q — 1 regular. 

• Then sequence of graphs G n , n — 2, 3, . . . converges to the infinite graph 
G = (V, E), where V — F x Z. The edges are either between the two 
vertices on the same level (U,k),k G Z, determined by D, or between 
the vertices of two consecutive levels (U, k) and (U, k + 1), given by the 
incidence matrix A in the way described above. 



ducible matrix. Let h G (p) be defined by (l2.1H2.2p . Then holds. 
(See Remark [5.21 ) 

In the example of G n = T((m', n)), n — 3, 4, . . discussed in the previous 
section, we have that U = T(m') and A is the identity matrix /. Hence A is 
also the identity matrix. 

7 The upper matching conjecture 

For r > 2 let K r ^ r be a complete bipartite graph on 2r vertices, where each 
vertex has degree r. Then 



P(t) : 



is the pressure of G. Assume that B(l) is an irre- 



(f)(l, K r>r ) = ( J H, / G Z +J and ip(x, K r , r ) = 





(7.1) 
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Conjecture 7.1 ( The upper matching conjecture.) Let G = (V,E) 
be a finite bipartite regular r -regular graph on 2qr vertices where 2 < q, r G N. 
Let qK rr be the graph consisting of q copies of K r r . Then 0(7, G) < <p(l, qK r ^ r ) 
for I — 0, . . . ,qr. 



In [TT] we proved the above conjecture for r = 2. We also showed that for 
r = 2 (f>(l, G) < 4>(l, qK 2 ^) for any 2 regular graph G on 4g vertices. {G does 
not have to be bipartite.) It is plausible that in the above conjecture one can 
drop the assumption that G is bipartite. For / = 0, 1 the above conjecture 
is trivial. For / = qr the above conjecture follows from the Mine conjecture 
proved by Bregman [3]. 

Let K(r) be an infinite countable union of K rr . Let hK( T ){p) be defined as 
in (l2.lH2.2p where G n = nK r ^ r , n G N. Let G n = (V n , E n ), n G N be a sequence 
of r regular bipartite graphs, where jfV n — > oo. Let h{c„}(p) be defined as 
in (I2.1H2.2I) . Assume for simplicity of the exposition that j^V n = 2g n r. Then 
Conjecture 17. II yields 0(/, G n ) < <p(l, q n K rjr ) for n G N. Hence the UMC yields 
the AUMC: h G {p) < h K{r) {p) for any pE [0,1]. 

We use the pressure PK(r)(f)i as pointed in Remark I5T2"1 to compute h K ^(p). 
Clearly the matching generating polynomial of qK r r is ip(x, K r ^ T ) q . Hence 

P K{r) {t)= gU= ll l} , teR. (7.2) 

This formula follows also from the results of the previous section, where F = 
K r ^ T and the incidence matrix A between two levels (U, 1) and (U, 2) is the 
zero matrix. Then p(B(t)A) = ip(e 2t , K r , r ). (}5T4"j) yields 

(p(t)) = ■cgSL.d)'^ _ (7 . 3 ) 



Vr n (:) 2 /!(2/)e 2Z * 
where pit) = ^ l=0 KlJ — V and * e R - 

2rELo(I) ^ 

Conjecture 7.2 ^ T/ie upper asymptotic matching conjecture.) Let 

Gk = (Vk,Ek),k £ N be a sequence of r regular bipartite graphs, where 
#Vfc — > oo. Lei /i{G fc }(p) defined as in $2.1W2.2\) . Let h K ^{p) be defined 
by (TJfy . Then h {Gk} (p) < h K{r) (p) for any p G [0, 1]. 

It is plausible to assume that Conjecture 17.21 holds under the assumption 
that each G n is an r-regular graph. 

Theorem 7.3 Let 2 < r G N and assume that G n = (V n ,E n ),n G iV is a 
sequence of r-regular bipartite graphs such that jfV n — > oo. Let h{c n }{p),P £ 
[0, 1] be defined by [2~NEEl . Then 

h{ Gn }(p) < min(upp pl (p), upp r]2 (p)) for all p G [0, 1], (7.4) 
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where 



plogrl 1 

u PPr,l(P) := — 2" 2 Pl ° gP_ ^ _P ) l0g ( 1 ^ 7 - 5 ) 

u PPr,2(p) : = ^ ^ r ~ ^(Pl0gp+ (1 -p)log(l -J))). (7.6) 

Proof. We claim that for any G = (V1UV2, E) G G(2m, r) the following 
inequality holds 

Indeed, let U\ C Vi be a subset of cardinality I. Consider the induced bipartite 
graph G(U\ U V2). Then G(J7i U Vj) is also an induced subgraph of the graph 
G' = (Vi U V 2 , E'), where the induced subgraph G'(Vi\Ui U V 2 ) is the complete 
bipartite graph on Vi\Ui U V 2 . It is straightforward to see that <f>(l, G') = (m — 
l)\<f>{l,G{Ui U V 2 ). The Bregman inequality [3] yields <fr(l, G') < (m!)^(r!)^. 
Since the number of choices of U\ is (™) we deduce the inequality (|7.7p . Let 
C n = (14, -E"n) G {?(2m n , r), n G N. Let /„, G [0, m n ]flZ be a sequence satisfying 
(I2.2p . Take the logarithm of the (17.71) divide by j^V n = 2fn n and let n — > 00 
to obtain that /i{G„}(^) < u PPr,i(p)- 
Our next inequality is 

0(Z,G)< (7)^ (7-8) 

Let C/i C Vi, = /. Since each vertex of G has degree r it follows that each 
vertex in Ui that (f>(l,G(Ui U Vi)) < r '- As one have (?) choices of the set U\ 
we obtain (17.81) . The above arguments imply that h{G n }{p) — u PPr2(p)- Hence 
(El holds. ' " □ 

Figure [3] gives the plot of h K ^, upp r 1 , upp r 2 for r = 4. From this graph and 
the corresponding graphs for r = 3, 8, which is not plotted here, we see that 
min(upp rl (p),upp rl (p)) — /ix(r)(p) decreases with r. Moreover the intersection 
point of the graphs upp r 1 and upp r 1 moves to the left as r increases. 



8 Computational results 

We have checked the asymptotic matching conjectures for several families like 
those described in Section [61 In each case we choose U to be a cycle Ci of length 
/ for several values of I and then varied the connection matrix A. In each case 
we used the described transfer matrix method to compute the entropy for 
several values of p and then compared with the conjectured bounds. In all 
cases the conjectures were found to hold. 

In order to test the conjectured lower bound for a given choice of U and A 
we first constructed the transfer matrix B(t) for the given graph. Given B(t) 
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we can directly compute P(t) from the maximum eigenvalue as in (15.21) . Next 
we computed P'(t), using the equality 



where rjf and i] 2 are the left and right eigenvectors of B(t), normalized by 
the condition r/J^ = 1. (This is a standard variational formula, e.g. [7].) 
From these values we now compute h G {p{t)) using (15. 4p . So for each value of 
t we get a pair (p(t),h G (p(t))) telling us the asymptotic pressure h G {p{t)) at 
the density p(t) = P'{t). To make all computations exact we chose e* to be 
rational numbers, which yielded rational values for all matrix entries. 

Example 8.1 (r=4) In our first family we let I, the length of the cycle 
U = C\, vary from 4 to 8. We tested all permutation matrices A, which give 
every vertex (u, k) in G n one neighbor in the level k — 1 and one in the level 
k + 1, and give rise to a bipartite G n . We thus have a family of bipartite 
4-regular graphs which includes the standard square lattice tori. 

In Figure^ we plot the difference between the actual values of the entropies 
ha{p), for all choices of A, and the lower asymptotic matching conjecture for 
a given range of densities p. The highest curve correspond to the normal torus 
graph, i. e it is the graph with the maximum number of matchings of a given 
size in this family. 

Example 8.2 (r=3) For our second example we again chose U = C\ to 
be a cycle of length I. Here we chose A so that if we number the vertices on 
the cycle the even vertices in the level k have an an odd vertex as a 

neighbor in the level k + 1 for k = 1, . . . , n. G n is cubic, and for I and n 
even G n is bipartite. This family includes the torus graphs obtained from the 
hexagonal lattice which all have girth at least 6. In this case we let the length 
of the cycle vary from A to 10 and again the conjecture was found to hold. In 
Figure O we display the difference between the actual values of the entropies 
h G (p), for all choices of A, and the lower asymptotic matching conjecture for 
a given range of densities p. Here the values typically stayed closer to the 
conjecture than for the 4-regular case, which is to be expected since this graph 
family has higher girth. 

Apart from the above tests, we also tested some more arbitrarily chosen 
connection matrices giving graphs of degree 6. This was done by U as a cycle 
and choosing the connection matrix A, having two l's in each row and column. 
Again the conjectures were found to hold but here the deviation up from the 
conjectured lower bound was even smaller. This is again expected since the 
conjecture should become more accurate for graphs of higher degree. 
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9 Infinite graphs with the maximal pressure 



In this section we give a partial justification for the computational result in 
Example 18.11 that the highest curve correspond to the normal torus graph. 

Theorem 9.1 Let F = (U,D) be an undirected graph and consider an 
infinite graph G = {V,E) defined as follows. V = U x Z, i.e. the vertices 
of G consists layers (U,k),k e Z. The edges E of G connect two vertices on 
the level k,j 6 Z only if \k — j\ < 1. The edges between any two vertices 
on level k are given by D. The edges between the level k and the level k + 1 
are given by #[/ x #[/ permutation matrix = {a uv ^) u ,v&u f or each k G Z. 
Thus ((u, k), (v, k + 1)) 6 E a uv ^ = 1. Let M(t), t e 1 and be the 

2# u x 2# u nonnegative matrices defined as in §5. Then the pressure Pa(t) of 
G is given as 

ij-+oo (i + j + l)#(7 

Lei Go 5e t/te infinite graph obtained by letting A^ to be the identity matrix for 
each k ZzTL. Then 

PcM = !25« > Po(() (9 . 2) 

/or any i el and any G o/ tae above form. In particular the monomer- dimer 
entropy h G of G, which is equal to -Pg(O), does not exceed ha = Pg o (0)- 

Proof. Fix i,j e N and let n = i + j + 1. Define G n = (14., i? n ) to 
be the following graph. V n consists of n copies of U, labelled as (U, k) for 
k = —j, . . . , i. The edges of E n are induced by the edges of G, except that 
the edges from the level i are connected to the edges of the level —j, which is 
identified with the level i + 1, by the connection matrix Aj. The arguments 
of §2 yield that i/j(e 2t ,G n ) = tr M{t)A^M{t)A_ j+1 . ..M(t)A t . Hence P G is 
given by fTI]) [13] ■ 

Consider now the case of G where A k = I. Then ( 19.11) yields -PgoW — 
l0S P #u • (^ ee ^ or exam pl e [9, §10] for the self-contained details of the ar- 
guments on matrices used here.) From the definition of M(t) it follows that 
M(t) is a nonnegative and symmetric matrix. Hence p(M(t)) = \\M(t)\\, 
where ||M(t)|| is the I2 operator norm of M(t). Since A& is a permutation it 
follows that At is also a permutation matrix. Hence \\Ak\\ = 1. Thus 

\\M(t)A_ j M(t)A_ j+1 ...M{t)A i \\ < \\M(t)\\ t+3+1 = p(M(t))^ +1 => 



tr M{t)A_jM{t)A_ j+1 . . . M(t)Ai < 2* u p{M{t)) i+j 



This proves fl9T2|) . 

From the definition of monomer-dimer entropy of G [12] it follows that 
h G = P G (0). Hence h G < h Go . □ 
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Conjecture 9.2 Let the assumptions of Theorem \9. II hold. Then for any 

pe[0,i] h G (p)<h Go ( P ). 
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r=4 
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Figure 3: hjc(4) -green, upp 4 1 -blue, upp 4 2 -orange, color online 
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Figure 4: Difference between actual entropy and the lower asymptotic match- 
ing conjecture for 4-regular graphs with U = G% 
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Figure 5: Difference between actual entropy and the lower asymptotic match- 
ing conjecture for 3- regular graphs with U = C w 
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